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Abstract 

One of the factors limiting electron mobility in supported graphene is remote phonon scattering. 
We formulate the theory of the coupling between graphene plasmon and substrate surface polar 
phonon (SPP) modes, and find that it leads to the formation of interfacial plasmon-phonon (IPP) 
modes, from which the phenomena of dynamic anti-screening and screening of remote phonons 
emerge. The remote phonon-limited mobilities for Si02, HfC>2, h-BN and AI2O3 substrates are 
computed using our theory. We find that h-BN yields the highest peak mobility, but in the prac- 
tically useful high-density range the mobility in Hf02-supported graphene is high, despite the fact 
that HfC>2 is a high-/i dielectric with low-frequency modes. Our theory predicts that the strong tem- 
perature dependence of the total mobility effectively vanishes at very high carrier concentrations. 
The effects of polycrystallinity on IPP scattering are also discussed. 
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I. INTRODUCTION 

Graphene, a single-layer of hexagonally arranged carbon atoms jlj], has been long consid- 
ered a promising candidate material for post-Si CMOS technology and other nano-electronic 
applications on account of its excellent electrical [2] and thermal transport [3] properties. 
In suspended single-layer graphene (SLG), the electron mobility has been demonstrated to 
be as high as 200,000 cm 2 V _1 s _1 J4|. However, in real applications such as a graphene field- 
effect transistor (GFET), the graphene is physically supported by an insulating dielectric 
substrate such as Si02, and the carrier mobility in such supported-graphene structures is 
about one order of magnitude lower J4|. This reduction in carrier mobility is further exac- 
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transport properties is a result of exposure to environmental perturbations such as scattering 
by charge traps, surface roughness, and remote optical phonons which are a kind of surface 
excitation. Such environmental effects are encountered in metal-oxide-semiconductor (MOS) 
structures 8|. Hess and Vogl first suggested that remote phonons [sometimes also known 
as Fuchs-Kliewer (FK) [9] surface optical (SO) phonons] can have a substantial effect on 
the mobility of Si inversion layer carriers [lOj. Fischetti and co-workers later studied the 
effects of remote phonon scattering in MOS structures and found that high-K oxide layers 
have a significant effect on carrier mobility in Si [8|| and Ge ll|. This method was later 
applied by Xiu to study remote phonon scattering in Si nanowires [l^j]. In Refs. [8] and 
121 ] it was found that the plasmons in the channel material (Si) hybridized with the surface 
polar phonons (SPP) in the nearby dielectric material to form interfacial plasmon-phonon 
(IPP) modes. This hybridization occurrence naturally leads to the screening/anti-screening 
of the SPP from the dielectric material. Scattering with these IPP modes results in a further 
reduced channel electron mobility in 2D Si and Si nanowires. 

Likewise in supported graphene, remote phonon scattering is one of the mechanisms 
believed to reduce the mobility of supported graphene, with the form of the scattering 
mechanism varying with the material properties of the dielectric substrate. Experimentally, 
hexagonal boron nitride (h-BN) has been found to be a promising dielectric material for 
graphene, and it is commonly believed that this is at least partially due to the fact that 



remote phonon scattering is weak with a h-BN substrate 



13jj. On other substrates such 
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as Si 
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Il5] and SiC 16|, LL7j, the mobility of supported graphene is lower. Thus, it is 
important to develop an accurate understanding of remote phonon scattering in order to 
find an optimal choice of substrate that will minimize the degradation of carrier mobility in 
supported graphene. 

Alth oug h the subject of remote phonon scattering in graphene 
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22j and carbon nan- 



20] ) has been broached in the recent past, the basic approach used in the afore- 



mentioned works does not deal adequately with the dynamic screening of the SPP modes. In 
graphene, dynamic screening of SPP modes has its origin in SPP-plasmon coupling, and the 
two time-dependent phenomena have to be treated within the same framework. Typically, 



the coupling phenomenon is ignored, and screenin g o 
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;he SPP modes is approximated with 



23] , which is adequate for the case of 



25] but can lead to a miscalculation of the scattering rates since the 



use of static screening underestimates the electron-phonon coupling strength [18|, especially 
for higher-frequency modes. The failure to incorporate correctly SPP-plasmon coupling into 
the approach means that the dispersion relation of the SPP (or, more accurately, of the IPP) 
modes is incorrect and that the dynamic screening of the remote phonons is not accounted 
for in a natural manner. 

To understand the screening phenomenon in our situation, let us first give a bird's eye view 
of the physical picture. This picture is somewhat different from what is found in the more 
familiar semiconductor-inversion-layer/high-K-dielectric geometry, since the absence of a gap 



in bulk graphene renders its dielectric response stronger and qualitatively different - a 
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metal-like, as testified by the presence of Kohn anomalies in the phonon spectra 
than the response of a two-dimensional electron gas. Graphene plasmons interact with the 
SPP modes through the time-dependent electric field generated by the latter, and the former 
are forced to oscillate at the frequency of the latter (oj). When u is less than the natural 
frequency of the plasmon (u p ) , i.e. u < u p , the electrons can respond to the SPP mode and 
screen its electric field. On the other hand, when u > u p , the motion of the plasmons lags 
that of the SPP mode, resulting in poor or no screening, or even in anti-screening, which can 
actually augment the scattering field [28]. In bulk SiC>2, the main TO-phonon frequencies 
are around 56 and 138 meV. At long wavelengths (A > 10 _8 m), the plasmon frequencies for 
a carrier density of 10 12 cm -2 are comparable or smaller than the TO-phonon frequencies. 
Thus, a TF-type approximation is inadequate especially for describing the screening (or 
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more accurately, the anti-screening) of the 138 meV TO phonon modes. Our calculations 
suggest that, contrary to what is found in the semi conductor /high- k case [8j and to the 
claims made in Ref. [181] . the higher-frequency SPP modes cannot be ignored despite their 
reduced Bose-Einstein occupation factors at room temperature. 

It is our intention in this paper to provide a systematic description of the coupling between 
the substrate SPP and the graphene plasmon modes, and relate this coupling to the dynamic 
screening phenomenon. Our theory can be generalized to graphene heterostructure such as 
double-gated graphene although this falls outside the scope of our paper and will be the 
subject of a future work. We begin by deriving our model of the IPP system. Its dispersion 
is then calculated from the model. The pure SO phonon and graphene plasmon branches are 
compared with the IPP branches. Also, we compute the electron-IPP and the electron-SPP 
coupling coefficients for different substrates (Si02, h-BN, Hf02 and AI2O3). We show that 
the IPP modes can be interpreted as dynamically screened SPP modes. Scattering rates are 
then calculated and used to compute the remote phonon-limited mobility/ipp for different 
substrates at room temperature (300 K) with varying carrier density. The temperature 
dependence of /ipp at low and high carrier densities is compared. Using the/ipp results, 
we analyze the suitability of the various dielectric materials for use as substrates or gate 
insulators in nanoelectronics applications. We also discuss the effects of polycrystallinity on 
remote phonon scattering. 

II. MODEL 

A. Coupling between substrate polar phonons and graphene plasmons 

Our approach to constructing the theoretical model of the coupled plasmon-phonon sys- 
tems follows closely that of Fischetti, Neumayer and Cartier {§} although some modifica- 
tions are needed to describe the plasmon-phonon coupling. One of the primary difficulties 
in describing the coupled system is the anisotropy in the dielectric response of graphene: 
graphene is polarizable in the plane but its out-of-plane response is presumably negligible. 
If the graphene sheet is modeled as a slab of finite thickness with a dynamic dielectric re- 
sponse in the in-plane direction [e| r (w) = ^gr{^ — u 2 /uj 2 ) where u p is the plasma frequency] 
and none in the out-of-plane direction \eg r {u) = constant], the dispersion of the SPPs re- 
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mains unchanged, indicating that the SPP and plasmon modes are uncoupled. This absence 



of couplin 
graphene 



; i s implicitly assumed in much of the current literature on SPP scattering in 
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29| although it has already been shown to be untrue in 2D Si [8] and 
Si uanowircs j ! 2|, Furthermore, there is considerable experimental support for the coupling 
of graphene plasmons to the SPPs 
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33] . As we will show later, accounting for this cou- 
pling results in the formation of IPP modes which are screened/anti-screened and scatter 
charge carriers in graphene more weakly/strongly than the unhybridized SPP modes. This 
'uncoupling problem' persists even when one inserts a vacuum region between the graphene 
slab and the substrate. Ultimately, this alleged lack of coupling can be traced back to the 
continuity of the electric displacement field D at the interface between the graphene slab 
and the substrate/vacuum. Given that the dynamic response of the graphene is only in 
the in-plane directions and that the coupling should be with the p-polarized waves of the 
substrate, the slab approach is not likely to be correct. To overcome this difficulty, we find 
it is necessary to treat the graphene as a polarizable charge sheet (as shown in Fig. [T]) rather 
than as a finite slab with a particular in-plane dielectric function. This polarization charge 
then generates a discontinuity in the electric displacement along the surface of the graphene. 
It is this discontinuity that couples the dielectric response of the substrate to that of the 
graphene sheet. The basic setup is shown in Fig. [TJ The graphene is an infinitely thin sheet 
co-planar with the x-y plane and floating at a height d above the substrate which occupies 
the semi-infinite region z < 0. Notation-wise, we try to follow Ref. 8|. In the direction 
perpendicular to the interface, the (ionic) dielectric response of the substrate is assumed to 



be due to two optical phonon modes, an approximation used in Ref. [8], that is: 

2 2 
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where Utoi and Uto2 are the first and second transverse optical (TO) angular frequencies 
(with utoi < UT02), and e^, e l ox and e® x are the optical, intermediate and static permittiv- 
ities. We can also express e ox (u) in the generalized Lyddane- Sachs- Teller form: 

p ,. a _oo Ko2-^ 2 )Koi 
tox{U) e 2 2 w 2 , ,2\ 

[UJ T02 - LO ) [iO TOl - UJ ) 

where u>loi and Ulo2 are the first and second longitudinal optical angular frequencies. The 
variables Q and R represent the two-dimensional wave and coordinate vector in the (x, y) 
plane of the interface, respectively. 




Figure 1: Schematic of set up of graphene-substrate system. The SLG is modeled as a infinitely 
thin (in the z-direction) layer of polarizable charge. A gap of d separates the graphene charge sheet 
and the substrate surface. 

As in Ref. we try to derive the longitudinal electric eigenmodes of the system since 
the transverse modes (given by poles of the total electric response) correspond to a vanish- 
ing electric field and so to a vanishing coupling with the graphene carriers. In effect, the 
longitudinal modes are the transverse-magnetic (TM) solutions of Maxwell's equations. It 
was also shown in Ref. js| that one may ignore the effects of retardation. Therefore, we need 
only to employ simpler electrostatics instead of the full Maxwell's equations. 

We begin our derivation by writing down the Poisson equation for the bare scalar potential 

<&, 

- V 2 $(R,z) = - Pox (R,z,t), (2) 

where p ox is the (periodic) polarization charge distribution at the surface of the substrate 
that is the source of scattering, and eo is the permittivity of vacuum. Equation (j2j) describes 
the electrostatic potential within the graphene. However, the effective scalar potential felt 
by the graphene carriers is different and should include the collective screening effect of the 
induced electrons/holes, which changes the RHS of Eq. (|2J). Hence, we modify Eq. (|2J) by 
adding a screening charge term on its RHS, and we obtain the Poisson equation for the 
screened scalar potential $ scr , 

- V 2 $ scr (R, z,t) = — [p ox (R, z, t) + p scr (R, z, t)) , (3) 
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where p scr is the screening charge term. The integral form of Eq. ([3]) is: 

$ scr (R, z, t) = $(R, z,t) + J dR'dz'G(Rz, RV)p scr (RV, t) (4) 

where G(R,z, RV) is the Green function that satisfies the boundary conditions [see 
Eqs. (Il4p ], and the equation: 

- V 2 [e(R, z)G(Rz, RV)] = 5(R - R', z — z') . (5) 

The bare potential $(R, z, t) is defined as $(R, z, t) = J dH'dz'G(Ti,z, R'z')p ox (R,'z' , t). The 
second term on the RHS of Eq. (j3J) represents the screening charge distribution. The bare 
and screened potentials can be written as sums of their Fourier components: 

<S>(R,z,t) = J2<pQAz)e- tQ - K e^ , (6a) 
Q 

$ scr (R, z,t) = Y, C^)e" lQ - R e^ (6b) 
Q 

where it must be understood that only the real part of Eq. (J6]) is to be taken here and in the 
following sections. Given the cylindrical symmetry of the problem, the Fourier components 
4>q iUJ and depend only on the magnitude of the wave vector Q. 

From Eq. (jlj) we obtain the following expression for the z-dependent part of the Fourier- 
transformed screened potential: 

C^K"* = thM** + / MG Q {z, z')p s ^(z', t) . (7) 

Equation (J7J) is solvable if the polarization charge pq^ is expressed as a function of the 
screened scalar potential. Here, we assume that pq^ responds linearly to (pq^, and write 
the screening charge term as: 

P&CM) = e 2 n(Q, W )/(z)0^(z)e^ (8) 

where H(Q, u) is the in-plane 2D polarization charge term, and f(z) governs the polarization 
charge distribution in the out-of-plane direction. For convenience, we model the graphene 
as an infinitely thin sheet of polarized charge and set f(z) = 5{z — d). Combining Eqs. (171) 
and (jSJ), we obtain the expression: 

4>%*( z ) = 4>qM + * [ dz'GQ&z'MQi^fiz'^iz') . (9) 
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The expression in Eq. ([9]) becomes: 



<VqM = <P Q ^z)+e 2 G Q {z,d)Tl{Q,u)<pQ:M 



+ e 2 G Q (z } d)n(Q, u^Q^d) + e A G Q {z, d)U(Q, u)G Q (d, d)U(Q, uj)<p s ^{d) 



, * v e 2 G Q (*,rf)n(Q,u;) 

and the corresponding component of the electric field perpendicular to the interface at z = 
is: 



z-E 



Q,LJ 



(z,t) dG Q (z,d) e 2 n(Q,w) 



G>2 



2=0 



92 l-e 2 G Q (d } d)Tl(Q } uj) 



4>q,M 



For notational simplicity, we write: 



VqM = 4>qM + G Q {z, d)V Q ^{d) 



2=0 

(10) 



'11a) 



z • E, 



,(z,t) dG Q (z,d) 



dz dz 
e 2 U(Q,cu) 



(lib) 



where 

Vq ' u = l-e^G Q (d,d)U(Q,u) ' 
Here, we emphasize that Eq. (11 lap is the key to determining the dispersion relation as we 



(12) 



shall show later. 

The Green function Gq(z,z') in Eq. (JTj) obeys the relation: 



-Q&" g2 ) G Q (z,z') = ±-5(z-z') 



(13) 



We require the Green function to satisfy the following conditions at and away from the 
interface {z — 0). 



dGniz = 0+ z') ^dGniz = 0-,z') 



e - 



dz 



dz 



(14a) 



G Q (z = + ,z f ) = G Q (z = 0-,z') 



(14b) 



G Q (z<0,d)=G Q (0,d)e +Qz 



(14c) 
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Gq(z > d,d) = G Q (d,d)e 



-Q(z-d) 



(14d) 



The solution to Eq. ( fT3l) is 



Gn(z, z) 



2e Q 



A _ Ae- Q l* +a/ l) ,z > 



^(l-A)e-^-'l 



z < 



(15) 



where 



A 



,oc 



The bare potential in Eq. (11 lap can be written as: 



Aie- Qz , 2 > 
A 2 e +Qz , z < 



where Ai and A 2 are the amplitudes of the bare potential for z > and z < respectively. 
Thus, the expression for the screened potential in Eq. (II lap is: 



A x e-^ z + C Q (z, rf)P Q , w A ie -« d , z > 
A 2e +^ + G Q (z, d)V Q ^A x e- Qd , z < 



(16) 



At the interface z = 0, the continuity of the component of the electric field parallel to the 
interface requires the continuity of (fiq^, i-e- = + ) = 4'Q r u( z = giving us: 



A x + A!G Q (z = 0+ d)V Q ^e~ Qd = A 2 + A x G Q {z = 0" djTW 



(17) 



Similarly, the continuity of the perpendicular component of the electric displacement, i.e., 

dz 



£o M^) =£OT ( lj) *S^) 1 e a d 8 to: 



dz 



<3 



dz 



1 dG Q (z = ,rf) Qd 



Q 



dz 



8) 



Substituting Eqs. (I14a[) and ()14b|) into Eqs. (ITT)) and (TT8)) . we obtain the following relations: 



Ax = A 2 



(19a) 



. l dG Q (z = 0+,d) Qd 

1 ~ ~ ' Q,ui^ 



Q 



dz 



l dG Q (z = 0-,d) Qd 
r „ „ / o.iiiC 



<9z 



(19b) 
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Rearranging the terms in Eq. (119b|) we obtain: 



eo + e ox (co) + [e ox (uj) - e~] G Q (0, d)V Q ^- Qd = 



(20) 



which can be rewritten as: 



[e + e m {u)\ {1 - [G Q {d, d) - G Q (0, d)e~ Qd ] e 2 IT(Q, w)}-(e + e£) G Q (0, d)e 2 U(Q, u)e~ Qd = 



or more explicitly: 

\e ox (u) + e ] 



2Qdj e 2 n(Q,o;) 



2e Q 



e 2 H(Q,cj) _ 2Qd 

Q 



o . 



(21) 



Equation (120]) gives us the dispersion of the coupled plasmon-phonon modes, and is some- 
times called the secular equation Physically, we expect three branches (two phonon and 
one plasmon). We write the coupled plasmon-phonon modes as 0Jq\ oJq^ and oj^q for each 
Q-point. In the limit d — > oo, Eq. ( J2TT) becomes: 

e 2 U(Q,u) 



[e ox (u) + e ] 



1 - 



2e Q 







which gives us as expected the dispersion for the two uncoupled SPP branches and the single 
plasmon branch in isolated graphene. 



B. Plasmon and phonon content 

The solutions of Eq. ( )2~Tj) represent excitations of the IPP modes. However, the effective 
scattering amplitude of a particular mode may not be substantial if it is plasmon-like. 
Scattering with a plasmon-like excitation does not necessarily lead to loss of momentum since 
the momentum is simply transferred to the constituent carriers of the plasmon excitation 
and there is no change in the total momentum of all the carriers. On the other hand, 
scattering with a phonon-like excitation does lead to a loss of momentum since phonons 
belong to a different set of degrees of freedom. Therefore, as in Ref. [8], it is necessary to 
define the phonon content [35] of each IPP mode. The phonon content quantifies the modal 
fraction that is phonon-like and modulates its scattering strength. Likewise, we can also 
define the plasmon content of the mode. To find the plasmon content, we first consider the 
two solutions wl 9 '°^ {a = 1, 2) obtained from the secular equation Eq. (l2TT) by ignoring the 
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polarization response [setting U(Q,uj) = 0]. Following Ref. [8j, the plasmon content of the 
IPP mode Wn is defined here as: 



(i)2 
CU>/ - UJ. 



(-9,1)2^ 



)(c« 2 -4-^ 2 ) 



("9,2)2^ 

Q 



f, ,«2 (i)2w (i)2 



(fe)2x 



where the indices k) are cyclical. Note that the expected 'sum rule' 



(22) 



(23) 



i=l 



holds. Equation ( |23|) implies that the total plasmon weight of the three solutions is equal to 
one (as it would be without hybridization). The (non-plasmon) phonon content is then 
defined as 1 — ^ 9 \uq > ). In order to distinguish the phonon-1 and phonon-2 parts of 
the non-plasmon content, we need to define the relative individual phonon content. For 
phonon-1, this is accomplished by ignoring its response and replacing e ox (uj) in Eq. (1211) 
with e^(w| 02 — uj 2 )/ {oJto2 ~ uj2 )- From the solutions of the modified secular equation 



(uJq TOl ' l > and w ^ T01 > 2 )^ relative phonon-1 content of mode i will be: 



( (i)2 (-T01,1)2 V# ,(i)2 



IK -w, 



(-TO!,2)2x 



(i)2 
Q 



, ,0')2\/ (*)2 w -5\ 



(24) 



where, as before, i, j and k are cyclical. The relative phonon-2 content can be similarly 
defined by replacing the superscript (— T01,a) with (— T02,a). Hence, the TO-phonon-1 
content will be: 



UJ, 



i?(™l) ( 4)) +J R(T02) (a; W ) 



1 



(25) 



The TO-phonon-2 content ^ T02S> (uj < q) can be similarly defined. Given Eqs. ( 1221) and (T55 
the following sum rules have been numerically verified: 

3 3 



E^ (roi) (^) = E* croa) ^g ) ) = 1 



i=l 



for each mode UJq 1 . 



$(9)( a; W) + $(TOl) i 



^ 02 \uj$) 
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C. Scattering strength 



As we have seen earlier, the IPP modes that result from the SPP-plasmon coupling have 
a different dispersion from that of the uncoupled SPP and plasmon modes. The electric field 
generated by the IPP modes is also different from that of the uncoupled SPP and plasmon 
modes. Since the remote phonon-electron coupling is derived from the quantization of the 
energy density of the electric field [8], we expect this difference in the electric field to be 
reflected in the scattering strength of the IPP modes. 

To find the scattering strength of an IPP mode, we have to determine the amplitude of its 
electric field. In Eq. (I16p there are three unknowns (A\, A 2 and u), two of which (Ai and u) 
can only be eliminated through Eqs. (1171) and (|l8p . To find Ai, we follow the semi-classical 
approach in Ref . [8j , where the time-averaged total energy of the scattering field is set equal 
to the zero-point energy. In the following discussion, we set A\ = Aq. We first compute the 
time-averaged electrostatic energy (UqQ associated with the screened field: 



(27) 



The angle brackets (. . .) denote time average. The volume integral in Eq. (127|) is the result 
of three contributions: one from the substrate (z < 0), one from the graphene-substrate 
gap (0 < z < d), and one from the region above the graphene (z > d). Each term can be 
converted into a surface integral. Adopting a 'piecewise approach' to evaluate the integral in 
Eq. (12 7J , we must evaluate three surface integrals. To do so, we need the explicit expressions 
for Gq(z, d): 

[1 - A) e +Q(- rf ) , z < 



i 



_i_ y e +Q{z-d) _ Ae -Q(,+d)j ; o < z < d 

_1_ [ e -Q(*-«0 _ \ e -Q(z+d)^ ? z > d 



(28) 



and for --^G Q (z, d): 



d 

-—G Q {z,d) 



_L ^ e +Q(z-d) _ Xe -Q(z+d)^ i0<z < d (29) 
J- [ e -Q(z-d) _ Ae -Q(*+d)] ? z > d 

We can now evaluate the electrostatic energy in the regions z<0,0<z<d and z > d. 



(U^ w ) = (U%,(z < 0)> + (U^(0 <z<d)) + (U%,(z > d)> 



(30) 
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As mentioned earlier, the volume integrals in Eq. ([27)1 can be recast as surface integrals. 
Thus, 

e QA 2 Q Q 



{U&{z < 0)> 



1 _ l dG Q (z = Q^,d) 



Q 



dz 



Q,ui e 



-Qd 



l + G Q (z = + ,d)V Q , UJ e- Qd 

(31a 



U^(0 <z<d)) 



1 1 0G Q (z = Q- d) . 



1 + G Q (z = (T, d):P Q)a ,e 



-Qd _ 1 = d + ,d) ^ Qd 







(9,2 



e -Q d + G Q (z = d + ,d)V Q ,ue 
(31b) 



' -Qd l dG Q (z=d ,d) Qd 



Q 



dz 



e- Qd +G Q (z = d-,d)V Q , UJ e- Qd 
(31c) 



In Eq. ( 131bj) . € ox (uj) is the dielectric function of the substrate, which we distinguish with 



the overhead bar, and distinct from e ox (u). As we shall see later, as in Ref. |8|, the function 
e ox (u) is chosen in a way consistent with the particular excitation that we want. Let us 
regroup the terms in Eqs. (131]) into those on the substrate surface at z = and those on the 
graphene at z = d. At z = 0, we have 
e VLA 2 o Q 



<w&(* = °)> 



2 

(4 ] )QAIQ 



while at z = d, we have: 



1 _ 1 dG Q (z = 0+,d) 



Q 



dz 



-Qd 



l + G Q (z = + ,d)V Q , U} e 



-Qd 



l dG Q (z = 0-,d) Qd 
+ Q dz 



l + G Q (z = 0-,d)V Q ^e- Qd 

(32) 



(U%,{z = d)) 



e QA 2 Q Q 
2 

e flA 2 Q Q 



Qd _ 1 dG Q {z = d + ,d) ^ 

1 ~ ~ ' QjU}^ 



Q 



dz 



-Qd 



Gq(z = d + ,d)V QiUJ e 



-Qd 



cm 1 dGn{z = d , <f) 



Q 



dz 



e- Qd + G Q (z = d-,d)V Q , UJ e- Qd 

(33) 



We have to be careful in computing (Uq 1 ^. Mathematically, it may seem that we ought to 
set (Uq r J) = (Uq^z = 0)") + (Uq 7 u (z = d)Y However, note that the term (Uq r J) accounts for 
the various excitation effects, ionic and electronic, but the term (UqJ^(z = d)) corresponds 
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to the charge singularity in the zero-thickness graphene sheet. It is a 'self-interaction' of the 
charge distribution in the graphene which has no dependence on uj, unlike (Uq^z = 0)), 
so that we should not expect it to contribute physically to the scattering of the graphene 
carriers. To see this more clearly, we rewrite Eq. 0271) as: 



= (lf dzdRD ■ e\ = /- J dzdR (e E 



where and P e are the polarization fields of the lattice (substrate) and the graphene 
electronic excitation respectively. The following identification can be made: 



= 0)> = (\ J dzdR (e E + P L ) ■ E 
(U%,(z = d)) = (± J dzdRP e -E 



and since we are only interested in the interaction of the lattice polarization field with the 
graphene carriers, we set \Un r J) = (Uq^(z = 0)), i.e.: 

= [eo + eoA^) + - C) G Q (z = 0, d)V Q ^ d ' 

x [l + G Q (z = 0,d)V Q , u e- Qd ] (34) 



where we have used the relationship ^^Gq(z = CT, d) = Gq(z = 0, d). In Eq. ( 134"1) . the first 
factor [eo + e ox (u) — . . .] resembles the secular equation, Eq. ( 1201) . Indeed, if we replace e ox (u) 
with e ox (co), then (Uq^ = as expected. This is no coincidence since UqI, represents the 
energy of the charge distribution present at the substrate- vacuum interface and the secular 
equation in Eq. (120|) is a statement about the absence of charges at the substrate-vacuum 
interface. This also confirms our earlier choice of excluding the contribution from the surface 
charges at z — d, since that contribution does not disappear when we replace e ox (u) with 
e ox (u). By regrouping the terms according to the position of their charge distribution in 
Eq. (130]) . we make the relationship to the secular equation Eq. (120]) manifest. Using Eq. (fT2]) . 
the expression for the screened electrostatic energy can be rewritten as: 

nA 2 Q Qv_ (i) w l (l-G Q (d,dyn(Q,u$) + G Q (0,dyn(Q,u$)e~^ 

X Q ' ' 2 L Q 1 V l-G Q (d,dyil(Q,u$) 

(35) 

We use the relationship (Wq^) = 2(Uq^) to obtain the time-averaged total energy, and 
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set it equal to the zero-point energy, i.e., \fvJJq = (Wq^,), so that: 



1 



1 - G Q (d, d)e 2 U{Q, u>$) + G Q (0, d)e 2 U{Q, u>$)e 



(i)\ p -Qd s 



l-G Q (d,d)e 2 U(Q,oo^) 



Therefore, the squared amplitude of the field is: 



A 2 = 



hoo 



(0 



1 - G Q (d, d)e 2 U{Q, oo^) + G Q (0, d)e 2 U{Q, u>$)e 



( i )\ p -Qd" 



2ttQ 



l-G Q (d,d)em(Q,uj$) 

To determine the strength of the scattering field, say for the TOl phonon, we take the 
difference between (i) the squared amplitude of the field with the TOl mode frozen and (ii) 
that of the field with the mode in full response. In (i), we set: 

/ 2 2 
-TOl.oo/ \ _ oo I U L02 ~ U 

^OX \~) S 



U! T02 — 00' 



and 



-TO 1,0, 



oo / U L02 



J LOl 



The squared amplitude of the TOl scattering field for oo = U0q is: 

A (O J*) 2 = ( 1 1 \ 

ZU( ^ \e x (uj ( q ') - e ox {uo y Q > ) e ox {uo y Q >) - e ox (uj K Q >) J 



$ (T01) (^) 



1 - G Q (d, d)e 2 U{Q, u>$) + G Q (0, d)e 2 U{Q, u%)e 



x 



(i)\ p -Qd s 



(36) 



l-G Q (d,d)e*IL(Q,u>$) 
The expression for the T02 scattering field can be similarly obtained. Therefore, the TOl 
effective scattering field can be written as: 

e 2 n(Q,c«) 



e- Qd + G Q (z,d) 



-Qd 



l-e 2 G Q (d,d)U(Q,uj§) 



(37) 



The scattering potential is: 

3 2 

V(R,z) = J2J2 eA ™i(Q^Q) 



1=1 n=l 



(38) 



where Qq (oq^) is the annihilation (creation) operator for the mode corresponding to Q 



and oo q . Generally, the graphene field operator can be written in the spinorial form as: 



*(*.*) = EE 



s=±l K 



se 



c K + 



s/C' 
C K 



e iK ' R v / 5(^ - d) (39) 
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where K.(K.') denotes the /C(/C) valley, and the +(— ) sign corresponds to the tt(tt*) band; 
c k( c k^) * s the annihilation (creation) operator of the s-band K electron state at the /C 
valley. Therefore, the interaction term is 



Hi 



hit 



dz J rfR^ f (R, z)V{R, z)tf(R, 



and, if we neglect the inter-valley terms, simplifies to: 



Hint ~ ^2^2 M Q a siK+Q,S2K (crJcjCk^ + SlS 2 C^ 



/C ' t c S2/c ' I a {l) + a (0t 

Q C K ) \ a Q + a -Q 



i=l si,s 2 K,Q 



(40) 



where 



^siKi,S2K2 



1 + sis 2 e"^ K i _eK 2) 



is the overlap integral that comes from the inner product of the spinors, and 



M. 



(0 



E 

u=l 



$ (TO ^t ) (w (0^ 



x 11 - G Q {d,d)e 2 IL{Q,u$) + G Q (0,d)e 2 U(Q,u^)e 



ox \ 

1 



1/2 



(41) 



is the electron-phonon coupling coefficient corresponding to the oo'q mode 



D. Landau damping 

At sufficiently short wavelengths, plasmons cease to be proper quasi-particle excitations 
because of Landau damping [281 ] . To model this phenomenon, albeit approximately, we 
take that to be the case when the pure graphene plasmon excitation, whose dispersion 
u = uj p (Q) is determined by the expression 1 — e 2 GQ(d, d)U(Q, u) = 0, enters the intra-band 



single-particle excitation (SPE) continuum 36j. This happens when the plasmon branch 



crosses the electron dispersion curve, i.e. when u p = hvpQ, and the wave vector at which 
this happens is Q c . When Q < Q c , the electron-phonon coupling coefficient in Eq. (140]) is 
that of Eq. (T4TT) . Although the lower-frequency IPP branches may undergo Landau damping 
from intra-band SPE as Uq 1 < VfQ, we still retain them because the sum rules in Eqs. ( 1231) 
and ( 1261) require us to maintain charge conservation {28 1. On the other hand, when Q > Q c , 
Landau damping is assumed to dominate all the IPP modes and the coupling between the 
substrate SPP modes and the graphene plasmons can be ignored. Instead of scattering with 
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three IPP modes for each given wave vector, we revert to using only two SPP modes. This 
allows us to satisfy the sum rules in Eq. (126]) . In this case, the electron-phonon coupling 
coefficient in Eq. (HIT) can be rewritten as: 



< = £ 

11=1 



e 2 hwX / 1 1 



2nQ y-TO»,oo {u V )+eo ,TO„0 (a; g) )+eo 

where / = SOi, SO2 indexes the SPP branch. 

III. RESULTS AND DISCUSSION 
A. Numerical evaluation 

Having set up the theoretical framework for electron-IPP interaction, we compute the dis- 
persion of the coupled interfacial plasmon-phonon modes and study the electrical transport 
properties. 

1. Interfacial plasmon-phonon dispersion 

In this section we compute the scattering rates from the remote phonons by employing 
the dispersion relation (oOq 1 ) and the electron-phonon coupling coefficient (Mq), which can 



be determined by solving Eqs. f[2"Ul) and PHI) , respectively. For simplicity, to solve the 
equations we use the zero-temperature, long-wavelength approximation for U(Q,u) 136 



atter 



37j: 



n(Q )W ) = ^ (42) 

where Ep is the Fermi level which can be determined from the carrier density n via the 
relation n = E F /(nh 2 v F ). 

In Fig. [2] we show the dispersion relation for an SiC>2 substrate with n = 10 12 cm -2 . The 
three coupled IPP branches are drawn with solid lines and labeled I, II and III while the 
dispersion of the uncoupled modes is drawn in dashed lines in the figure. The branches 
labeled 'SOV (61 meV) and 'SOV (149 meV) have a flat dispersion and are determined from 
the quantity: 

e + e ox (u) = 
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while the branch labeled 'Pure plasmon' is determined from the zeros of the equation: 

l-G Q (d,d)e 2 U(Q,cu) (43) 



which gives the dispersion of the pure graphene plasmons when the frequency dependence 
of the substrate dielectric function is neglected and only the effect of the substrate image 
charges is taken into account. We observe that in the long wavelength limit (Q — > 0), 
branches I, II and III converge asymptotically to the 'pure' plasmon, SOi, and SO2 branches 
respectively. On the other end, as Q — > 00, branches I, II and III converge asymptotically to 
the pure SOj, SO2 and plasmon branches respectively. At intermediate values of Q the IPP 
branches are a mixture of the pure branches. The coupling between pure SO phonons and 
graphene plasmons has often bee n ig nored in transport studies based on the dispersionless 



unscreened, decoupled SO modes 



18 



19 



21-23 



techniques, Hwang, Sensarma and Das Sarma 



29 



38 



On the other hand, using many-body 



39] have studied the remote phonon-plasmon 



coupling in supported graphene and were able to reproduce the coupled plasmon-phonon 



dispersion observed by Liu and Willis 



30 



3l| in their angle-resolved electron-energy-loss 



spectroscopy experiments on epitaxial graphene grown on SiC. Similar results of strongly 



coupled plasmon-phonon modes were reported by Koch, Seyller and Schaefer 



Fei and 



co-workers also found evidence of this plasmon-phonon coupling in the graphene-Si02 system 
in their infrared nanoscopy experiments 



32|. Given the increasing experimental support for 
the hybridization of the SPPs with the graphene plasmons, it is interesting to investigate 
the effect of these coupled modes on carrier transport in graphene. 



2. Electron-phonon coupling 

Here, the electron-phonon coupling coefficients Mq 1 of the IPP and the SPP modes are 
compared. Recall that IPP modes are formed through the hybridization of the SPP and 
graphene plasmon modes, and their couplin g to the graphene electrons are different to that 



of the SPP modes. It is sometimes assumed 
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2l| that the SPP modes are screened by the 



plasmons, and the IPP-electron coupling is weaker than the SPP-electron coupling. As we 
have discussed above, this assumption does not hold when the frequency of the IPP mode 
is higher than the plasmon frequency. We plot the Mq'Q for the SPP and IPP modes in 
Fig. |3jWe first notice that at small Q, the coupling terms for branches I and II are actually 
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Figure 2: Dispersion relation of coupled interfacial plasmon-phonon system with n = 10 cm . 
The three hybrid IPP branches are labeled I (blue), II (green) and III (red). The uncoupled pure 
SO phonon (1 and 2) and plasmon branches are drawn with dashed lines. In the limits Q — »■ and 
Q — > oo, the IPP branches converge to the pure phonon and plasmon branches. In between, they 
are a mix of the pure branches. 

larger than those for SOi and SO2, even though I and II are phonon-like. This is because at 
long wavelengths, u p < Uq for I — 1, 2, resulting in anti-screening, effect which enhances the 
SPP electric field. For the plasmon-like branch III, Mq 11 ^ is actually much larger than the 
those for SOi and SO2 over the entire range of Q values. When we take Landau damping 
into account, we use the coupling coefficients (shaded in gray in Fig. [3]) of branches I, II and 
II for Q < Q c and of SOi and SO2 for Q > Q c . 
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Figure 3: Plot of MqQ for n = 10 12 cm~ 2 in SiC>2- The IPP branches are labeled I (blue), II 
(green) and III (red), and the SPP branches are labeled SOi and SC>2- The cutoff wave vector Q c 
is drawn in gray dashed lines. When Q < Q c , we use the part of the IPP branches shaded in gray, 
and when Q > Q c , we use the part of the SPP branches shaded in gray. 

B. Substrate-limited mobility 

The momentum relaxation rate for an electron in band s with wave vector K can be 
written as: 



where N B (co) = ( e ^/ fesT - l)" 1 , f(E) = [ e (E~E F )/k B T + and E ^ = s ^|k|. In assuming 



the latter expression, we use the Dirac-conical approximation. Equation flU]) automatically 
includes the Fermi-Dirac distribution of the final states and remains applicable when the 
doping level is high. The individual scattering rates for the screened (I, II and III) and 
unscreened (SOi and SO2) branches at the carrier concentration of n — 10 12 cm -2 in SiC>2 
and Hf02 are plotted in Fig. HI Landau damping is taken into account by setting the coupling 



s, K) = — ^2 ^2 \ M Q a ^+Q,s'K [1 - ss' cos(6»k+q - 0k)] 



I s' Q 




(44) 
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Figure 4: Plot of scattering rates at n = 10 12 cm' 2 for different substrates: (left) Si02 and (right) 
HfC>2- In SiC>2, the plasmon-like branch III dominates the scattering rate at E = Ep. In HfC>2, 
branches III, SOi and SO2 dominate the scattering rate at E = Ep = 117 meV. The SPP branches 
(SOi and SO2) do not contribute much to the Fermi-level scattering rate in SiC>2 because of their 
higher frequencies and smaller occupation factors. 

coefficient of the IPP (SPP) modes to zero when Q < Q c (Q > Q c ). We observe that at low 
energies, the IPP scattering rates are much higher than the SPP ones. At higher energies, 
the SPP scattering rates increase rapidly. The dominant scattering mechanism around the 
Fermi level appears to be due to the plasmon-like branch III in SiC>2 and HfC>2. In addition, 
at the Fermi level in Hf02, the SPP branches have scattering rates comparable to those of 
branch III. This explains why the low density mobility of Hf02 is less than that of SiC>2. 
The expression for the IPP/SPP-limited part of the electrical conductivity is: 



°RP = ^f kBT J o f(E - E F )[1 - f(E - E F )]T tr (E)~ 1 EdE . (45) 

where g s = 2 and g v = 2 are the spin and valley degeneracies respectively. Only the 
contribution from the conduction band is included in Eq. (I45jl . We use Eqs. (I44j) and (1431) 
to compute the IPP/SPP-limited electrical conductivity by setting: 



RP 



>,K). 



(46) 



In making this approximation, we ignore the other effects (ripples, charged impurity, acoustic 
phonons, optical phonons, etc). The scattering rates from the acoustic and optical phonons 
tend to be significantly smaller and are not the limiting factor in electrical transport in 
supported graphene (40|. Impurity scattering tends to be the dominant limiting factor, but 
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its effects can be reduced by varying fabrication conditions. Thus, the conductivity using 
Eq. (146]) gives us its upper bound. We calculate the remote phonon-limited mobility as: 

aRP S A1\ 

mp = — (47) 

en 



where n = J °° f(E — Ep)EdE is the carrier density. For n = 10 12 cm -2 in SiC>2, we 

obtain fi^p w 40,000 cm 2 V _1 s _1 . This is more than the corresponding values reported in 
the literature (~ 1000 — 20, 000 cm 2 V -1 s -1 ) \^ although we have to bear in mind that 
it is an upper limit. Nonetheless, it suggests that IPP/SPP scattering imposes a bound on 
the electron mobility. 



C. Mobility results 



Although suspended graphene has an intrinsic mobility limit of 200, 000 cm V s at 



room temperature 



Pi 



typical numbers for graphene on Si02 tend to fall in the range 1000- 



20,000 cm 2 V x s 1 [14]. One significant reason for this drastic reduction in mobility is be- 
lieved to be the presence of charged impurities in the substrate which causes long-range 



Coulombic scattering 
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42fl and much effort has been directed towards the amelio- 



ration of the effects of these charged impurities. For example, it has been suggested that 
modifying the dielectric environment of the graphene, either through immersion in a high-K 
liquid or an overlayer of high-K dielectric material, can lead to a weakening of the Coulombic 
interaction and an increase in electron mobility 42J. On the other hand, actual experimental 



evidence in favor of this theory is ambiguous. Electrical conductivity data from Jang and 
co-workers 42| as well as Ponomarenko and co-workers 43] indicate a smaller-than-expected 



increase in mobility when a liquid overlayer is used. This suggests that mechanisms other 
than long and short-range impurity scattering are at play here. Here, we turn to the problem 
of scattering by IPP modes. 



1. Comparing different substrates 

Having set up the theoretical framework in the earlier sections, we now apply it to the 
study of the remote phonon-limited mobility of four commonly-used substrates: Si02, Hf02, 
h-BN and AI2O3. Silicon dioxide is the most common substrate material while Hf02 and 
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AI2O3 are high-K dielectrics commonly used as top gate oxides |38l . |44| . Hexagonal boron 
nitride shows much promise as both a substrate and a top gate dielectric material The 
study of the remote phonon-limited mobility in these substrates allows us to understand 
how electronic transport in supported graphene depends on the frequencies and relative 
permittivities of the substrate phonons. 

From Eq. (T4"5l) with the effects of Landau damping taken into account, we compute 
the remote phonon-limited mobility numerically, using the well-known Gilat-Raubenheimer 
method 45] to discretize the sum in Eq. (j45l) . We plot /i^p as a function of carrier density 
(n = 0.3 x 10 12 cm" 2 to 5.2 x 10 12 cm" 2 ) at 300 K in Fig. EJ Note that the mobility values 
of the high-K substrates (Hf02 and Al 2 03) are substantially lower compared to SiC>2 and 
h-BN in the carrier density range n < 2.0 x 10 12 cm -2 . Similar results have been found in 
MOS systems Jf|. Hexagonal BN has the highest mobility at low carrier densities because 
of its high phonon frequencies, which corresponds to low Bose-Einstein occupancy, as well 
as its weak dipole coupling to graphene. In general, fipp for all four substrates increases 
with n because the dynamic screening effect becomes stronger at higher carrier densities. 
At low carrier densities, the mobility is low for all the substrates because there is a large 
proportion of plasmons modes whose frequencies are lower than the SPP mode frequencies. 
Thus, their coupling to the SPP modes results in the formation of anti- screened IPP modes 
that couple more strongly to the carriers, a phenomenon that has been studied for polar 
semiconductors [28]. However, as n increases, the mobility for all four substrates rises 
because the plasmon frequency scales as co p oc n 1//4 , resulting in higher- frequency plasmon 
modes. Thus, the plasmon-phonon coupling forms screened IPP modes that are weakly 
coupled to the carriers. Furthermore, at higher carrier densities, Landau damping becomes 
less important as a result of the increasing magnitude of the plasmon wave vector Q c . 
Contrary to expectation, we find that the mobility for Hf02 exceeds those of other substrates 
at larger densities (n = 5x 10 12 cm" 2 ). At n = 5 x 10 12 cm -2 , Hf02 has the highest remote- 
phonon mobility followed by h-BN, Si02 and AI2O3. This is because the proportion of 
screened IPP modes increases with increasing carrier density. Given the small values of 
ojtoi an d u)to2 f° r HfO"2, its coupling coefficients are smaller as a result of stronger dynamic 
screening. This weaker coupling compensates in part the higher occupation factors. In 
contrast, the larger values of ojtoi and UT02 for h-BN imply that screening does not play 
a significant role at low carrier densities. Hence, its coupling to the graphene carriers does 
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not diminish as rapidly as carrier density increases. The computed (jlrp values for H1O2 and 
h-BN highlight the role of low-frequency excitations in carrier scattering. The low-frequency 
modes are highly occupied at room temperature and induce carrier significant scattering at 
low n. At higher n when dynamic screening becomes important, the low-frequency modes are 
more strongly screened and their coupling to the carriers becomes diminished more rapidly 
than that of high-frequency modes. 



2. Dynamic screening effects 

To compute the mobility for the case without any screening or anti-screening effects, the 
Landau damping cutoff wave vector is decreased, i.e., Q c — > 0, resulting in the replacement of 
all the IPP modes with SPP modes. We plot the SPP-limited mobility as a function of carrier 
density in Fig. [5] (solid symbols), and compare these results for the IPP-limited mobility. The 
SPP-limited mobility for different substrates spans a range of values varying over nearly two 
orders of magnitude. In the absence of dynamic screening or anti-screening, the SPP-limited 
mobility for HfC>2 is only around 1000 cm 2 V _1 s _1 at n = 10 12 cm -2 , more than an order of 
magnitude smaller than the corresponding IPP-limited mobility, because of its low phonon 
frequencies. This result is also clearly inconsistent with experimental observations, since 



significantly higher mobility values have been reported for HfCVcovered graphene |38l |46|. 
The drastic reduction of the computed mobility suggests that screening is very important 
for the determination of scattering rates in a coupled plasmon-phonon system with low 
frequency modes. In contrast, h-BN gives an SPP-limited mobility of ~ 110, 000 cm 2 V _1 s _1 
at n = 0.3 x 10 12 cm~ 2 , which is still close to the IPP-limited mobility, indicating that its 
high frequency modes are relatively unaffected by screening. The maximum SPP-limited 
mobility for AI2O3 is around 8, 400 cm 2 V _1 s _1 at n = 0.3 x 10 12 cm~ 2 , which is much smaller 
than the 19000 cm 2 V _1 s _1 extracted by Jandhyala and co-workers |47fl who used AI2O3 for 
their top gate dielectric. This disagreement reinforces the necessity of including dynamic 
screening effects. Furthermore, the carrier density dependence of SPP-limited mobility is 
different from that of IPP-limited mobility. The IPP-limited fipp increases rapidly with 
carrier density because dynamic screening becomes stronger at higher n, an effect that is 
not found in SPP-limited mobility. In contrast, SPP-limited fipp decreases monotonically 
with increasing n. 
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e ox ( e o) 


3.90 


5.09 


22.00 


12.35 
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e S ( e o) 


2.50 
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5.03 
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utoi (meV) 


55.60 


97.40 


12.40 


48.18 


w T 02 (meV) 


138.10 


187.90 


48.35 


71.41 



Table I: Parameters [see Eq. ([TJ] used in computing dispersion relation and scattering rates for 
Si0 2 , h-BN, Hf0 2 and A1 2 3 . They are taken from Refs. [8] and [29 1. 



Our results suggest that Hf0 2 remains a promising candidate material for integration with 



2l| 



graphene since its high static permittivity can reduce the effect of charged impurities 
while its IPP scattering rates are relatively low when the carrier density is high. Although 
its surface excitations are low-frequency, which results in high Bose-Einstein occupancy, 
this is offset by its relatively strong dynamic screening at higher carrier densities. Thus, 
IPP scattering does not represent a problem for its integration with graphene field-effect 
transistors. As expected, h-BN is also a good dielectric material since its high phonon 
frequencies imply a low Bose-Einstein occupation factor. Furthermore, its smooth interface 
results in a smaller interface charge density and is less likely to induce mobility-limiting 
ripples in graphene. 



3. Temperature dependence 

Remote phonon scattering exhibits a strong temperature dependence - stronger than for 
ionized impurity scattering - because the Bose-Einstein occupation of the remote phonons 
decreases with lower temperature. This change in the distribution of the remote phonons 
(IPP or SPP) necessarily implies that the electronic transport character of the SLG must 
change with temperature. At lower temperatures, scattering with the remote phonons de- 
creases, resulting in a higher remote phonon-limited electrical mobility. The dependence of 
the change in mobility with temperature is related to the dispersion of the remote phonons 
and their coupling to the graphene electrons. By measuring the dependence of the mobility 
or conductivity with respect to temperature, it is possible to determine the dominant scatter- 
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Figure 5: Calculated conductivity remote phonon-limited mobility for different values of carrier 
density and different substrates (Si02, Hf02, h-BN and AI2O3) at room temperature (300 K). The 
IPP-limited mobility values are plotted using solid lines with unfilled symbols while the SPP-limited 
mobility values are plotted using dotted lines with solid symbols. 

ing mechanisms in the supported graphene. Given that our model of electron-IPP scattering 
differs from the more common electron-SPP scattering model, comparing the temperature 
dependence of the substrate-limited mobility can enable us to distinguish between the two 
models. 

The mobility of supported graphene over the temperature range of 100 to 500 K for the 
different substrates is computed at carrier densities of n = 10 12 cm -2 and n = 10 13 cm~ 2 . 
For the purpose of comparison, we perform the calculation for the case with screening (IPP) 
and without screening (SPP). The results (1///rf vs. T) are shown in Fig. [6^ and b. In 
Fig. Ek, we plot the IPP- and SPP-limited inverse mobility at n — 10 12 cm -2 . As expected, 
the substrate-limited mobility decreases with rising temperatures for both the screened and 
unscreened cases. From the plots, we observe that there exists an 'activation' temperature 
for each substrate at which the inverse mobility increases precipitously. For SiC>2, that 
temperature is around 200 K in the screened case and around 120 K in the unscreened case. 
This difference is striking and may be used to distinguish the IPP model from the SPP model 
at low carrier densities. In all four substrates, the slope of 1/hrp with respect to T is also 
steeper in the IPP-limited case than in the SPP-limited case. In Fig. |6b, we plot again the 
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IPP- and SPP-limited inverse mobility but at a much higher carrier density of n = 10 13 cm -2 . 
The IPP-limited 1//xrp is about three orders of magnitude smaller than the SPP-limited 
1/Hrp from 100 to 500 K. At high carrier densities, IPP scattering is insignificant and any 
changes in total mobility with respect to temperature cannot attribute to IPP scattering. 

The results in Fig. suggest that if IPP modes are the surface excitations that limit carrier 
transport in SiCVsupported graphene at room temperature, then the mobility would have 
a significant increase at around 200 K for n = 10 12 cm -2 . However, this IPP temperature 
dependence disappears at much higher carrier densities (~ n = 10 13 cm -2 ) because the IPP 
coupling to electrons becomes so weak that it no longer contributes significantly to carrier 
scattering. The results in Fig. also shows that the [Xrp increases monotonically with 



n. This should be contrasted with the result of Fratini and Guinea [18| who found that 
IIrp decreases as ~ a ^ room temperature. This is because the coupling coefficient 

limg^oo Mq, which is proportional to the matrix element, scales as \j\J~Q in the SPP model 
with static screening. In Fig. [3j linig^oo MqQ decreases with Q, implying that liniQ^oo Mq 
scales as Q a where a < — 1. In other words, the coupling coefficient vanishes more rapidly 
with Q in the IPP model than the SPP model. Our fipp results parallel those in Ref. 



48| in which the remote phonon-limited mobility increases with the carrier mobility in a 



2-dimensional electron gas system in the Si inversion layer with high-K insulators. 

In supported SLG, the carrier mobility is limited by three scattering mechanisms: long- 
range charged impurity, short-range defect and remote phonon scattering J49 1. The intrinsic 
phonon scattering processes in graphene can be effectively neglected. Of the three scattering 
mechanisms, only remote phonon scattering is strongly temperature dependent. The IPP 
model suggests that remote phonon scattering diminishes with increasing carrier density. 
Thus, the experimental consequence is that the temperature dependence of the mobility 
in supported-SLG should weaken at higher carrier densities. On the other hand, the SPP 
model predicts that the temperature dependence of the mobility should increase at higher 



carrier densities 



18]. This difference in the temperature dependence of the total mobility 



between the two models should be easily discriminable in experiments. 
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Temperature 7~(K) 

Figure 6: Inverse SPP and IPP-limited mobility versus temperature at (a) n = 10 12 cm -2 and (b) 
n = 10 13 cm -2 for SiC>2, Hf02, h-BN and AI2O3. is strongly temperature dependent at low 

carrier densities only. 

4- Disordered graphene 

We discuss qualitatively the interfacial plasmon-phonon phenomenon in disordered 
graphene. It is well-known that graphene grown by chemical vapor deposition (CVD) j^] 
is generally polycrystalline and contains a high density of defects. In supported graphene, 
charged impurities from the substrate and other defects scatter graphene carriers. These 
defects can affect the dynamics of plasmons in graphene which may in turn affect the hy- 
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bridization between the plasmon and the SPP modes. At short wavelengths, the plasmon 
lifetime rapidly decreases as a result of Landau damping which results in the decay of the 
plasmons into single-particle excitations. At long wavelengths, the plasmon lifetime can 
be affected by defects in the graphene. As far as we know, there is no theory of graphene 
plasmon damping from defects. However, it has been pointed out that long-wavelength plas- 
mons in polycrystalline metal undergo anomalously large damping due to scattering with 



structural defects 



51] . If this is also true in polycrystalline or defective graphene, then it 



implies that the long-wavelength surface excitation in supported graphene are SPP, not IPP, 
modes. 

To model phenomenologically this damping of long-wavelength plasmon modes in poly- 
crystalline graphene with defects, we set another cutoff wave vector Qd below which the 
surface excitations are SPP and not IPP modes. Qd is possibly related to the length scale A 
of the inhomogeneities or defects in graphene. As a guess, we choose A= 6 nm, which is a typ- 
ical autocorrelation length of 'puddles' in neutral supported graphene 52], and set Qd = 1/A. 
Hence, in our model, for Q > Q c and Q < Qd, the surface excitations are SPP modes while 
for Qd < Q < Qc, they are IPP modes. We compute the remote phonon-limited mobility at 
300 K and plot the results in Fig. [7J We find that the long-wavelength SPP dramatically 
alters the carrier dependence of Hrp in SiC>2, HfC>2 and AI2O3. In perfect monocrystalline 
graphene, fiptp reaches ~ 2 x 10 6 cm 2 V _1 s _1 in HfC>2 and SiC>2 at n = 5 x 10 12 cm~ 2 . On 
the other hand, in polycrystalline graphene with defects, it drops to the range of 10 4 to 
10 5 cm 2 V -1 s -1 . For h-BN, /i#p is quite relatively unaffected by the long- wavelength SPP 
modes except at low carrier densities (n < 0.5 x 10 6 cmT 2 ). 

This change in remote phonon-limited mobility highlights the possible role of defects 
in the surface excitations of supported graphene. We emphasize that our treatment is 
purely phenomenological and a more rigorous treatment of plasmon damping is needed 
in order to obtain a more quantitatively accurate model. Nevertheless, it emphasizes the 
relationship between dynamic screening and plasmons. In highly defective graphene, the 
surface excitations may be unscreened SPPs rather than IPPs because of plasmon damping. 
This should be taken into account when interpreting electronic transport experimental data 
of exfoliated and CVD- grown graphene. 
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Figure 7: Remote phonon-limited mobility in perfect monocrystalline (clear symbols) and defective 
polycrystalline graphene (solid symbols) for Si02, HfC^, h-BN and AI2O3. As carrier density 
increases, [Irp also increases. The use of long-wavelength SPP modes leads to a significant decrease 
in hrp in polycrystalline graphene. 

IV. CONCLUSION 

We have studied coupled interfacial plasmon-phonon excitations in supported graphene. 
The coupling between the pure graphene plasmon and the surface polar phonon modes of 
the substrates results in the formation of the IPP modes, and this coupling is responsible for 
the screening and anti-screening of the IPP modes. Accounting for these modes, we calculate 
the room temperature scattering rates and substate-limited mobility for SiC>2, HfC>2, h-BN 
and AI2O3 at different carrier densities. The results suggest that, despite being a high-K 
oxide with low frequency modes, HfC>2 exhibits a substrate-limited mobility comparable 
to that of h-BN at high carrier densities. We attribute this to the dynamic screening of 
the HfC>2 low-frequency modes. The disadvantage of the higher Bose-Einstein occupation 
of these low-frequency modes is offset by the stronger dynamic screening which suppresses 
the electron-IPP coupling. Our study also indicates that the contribution to scattering 
by high-frequency substrate phonon modes cannot be neglected because of they are less 
weakly screened by the graphene plasmons. The temperature dependence of the remote 
phonon-limited mobility is also calculated within out theory. Its change with temperature is 
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different at low and high carrier densities. We find that in the IPP model, the temperature 
dependence of the mobility diminishes with increasing carrier density only, in direct contrast 
to the predictions of the more commonly used SPP models. The implications of the damping 
of long- wavelength plasmons have also been studied. We find that the it leads to a substantial 
reduction in the remote phonon-limited mobility in SiC>2, HTO2 and AI2O3. 
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